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Abstract 

An exact analogue of the method of averaging in classical mechanics is con- 
structed for self-adjoint operators. It is shown to be completely equivalent to the 
usual Rayleigh-Schrodinger perturbation theory but gives the sums over interme- 
diate states in closed form expressions. The anharmonic oscillator and the Henon- 
Heiles system are treated as examples to illustrate the quantum averaging method. 

PACS Code: 03.65.-W, 31.15.+q, 02.30.Mv, 02.90.+p 

1 Introduction and motivation 

The failure to obtain exact solutions for most mechanical systems of interest (e. g. plan- 
etary motion) has prompted the search for perturbation techniques almost immediately 
after the conception of Newtonian mechanics (see for some history on the subject). 
At about the turn of the last century Lindstedt 0, Poincare 0, and later von Zeipel 
developed a perturbation method for classical Hamiltonian systems using an averaging 
procedure in phase space. Despite its lack of convergence in many cases this method 
which we shall henceforth refer to as the Poincare-von Zeipel method has been a widely 
used one since it yields at least asymptotic expansions. 
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Concerning the ability to find exact solutions nothing much changed with the advent of 
quantum mechanics. There it turned out to be equally important to develop perturbation 
methods and this was done simultaneously with the beginning of quantum mechanics 
by Schrodinger [Q. Due to previous contributions to similar perturbation techniques in 
other wave equations by Lord Rayleigh this theory has been named Rayleigh-Schrodinger 
perturbation theory and has later been given a rigorous mathematical basis in the work 
of Kato i and Relhch §. 

In this paper it is shown that the two methods are identical. More precisely, it will be 
shown that in quantum mechanics an exact analogue of the classical Poincare-von Zeipel 
method can be formulated with the help of an averaging techique for self adjoint operators 
analogous to the classical method and that the resulting quantum Poincare-von Zeipel 
perturbation theory is identical to the Rayleigh-Schrodinger theory. 

The analogy between the cassical and quantum case is based entirely on the structure 
of the equations appearing in the algorithm and the strucure of the method used to 
solve them (averaging). The starting point in the classical case is a Hamiltonian function 
on phase space with a perturbing Hamiltonian, whereas in the quantum case we start 
from a self adjoint Hamiltonian operator with a perturbing operator. If and how the 
two Hamiltonians are related is of no interest here. No quantization or other quantum- 
classical map (e. g. semi-classical correspondence) is needed to construct the quantum 
analogue of the classical Poincare-von Zeipel theory. 

Viewing Rayleigh-Schroding quantum version of the classical Poincare-von 

Zeipel method yields (apart from a purely conceptual viewpoint) one possible advantage: 
it gives closed expressions for the sums over intermediate states which appear in the 
corrective terms for the eigenvalues and eigenfunctions. 

The method of quantum averaging has also been used to construct a quantum analogue 
of Kolmogorov's superconvergent perturbation theory This new quantum "supercon- 
vergent" perturbation theory is substantially different and from the standard Rayleigh- 
Schrodinger method and initial numerical studies in some examples indicate much better 
convergence properties 0. While using quantum averaging to construct analogues of the 
classical Poincare-von Zeipel and superconvegent methods they needed to be compared 
with existing perturbation methods in quantum mechanics. In |^ we have shown that the 
quantum superconvergent method yields a new kind of perturbation theory and in this 
paper we show that the quantum Poincare-von Zeipel method is identical to the standard 
Rayleigh-Schrodinger theory. 

In classical mechanics the Poincare-von Zeipel series is in most cases divergent and 
yields only an asymptotic series. This is similar in its quantum equivalent the Rayleigh- 
Schrodinger series. Since we prove the equivalence of the two methods we shall not state 
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all technical details necessary to make all steps rigorous but refer the reader to the vast 
mathematical literature dealing with the Rayleigh-Schrodinger (or Kato-Rellich) pertur- 
bation theory (see e. g. fT^)- Just as certain quantities diverge in standard Rayleigh- 
Schrodinger theory (in cases where convergence conditions fail) the power expansions in 
e formally written down here may not converge in which case the perturbation algorithm 
gives only asymptotic information, all sums have to be replaced by finite ones up to N, 
and equations have to be read modulo 0(e^"^^) for any finite A^. 
The paper is organized as follows: 

In section ^ we present the classical Poincare-von Zeipel perturbation theory and 
method of averaging in such a way that it can easily be generalized to quantum mechnics 
which is done in section |^. 

In section ^ we apply the quantum Poincare-von Zeipel and averaging method to a 
Hamiltonian with pure point spectrum, show that up to second order all results from the 
Rayleigh-Schrodinger theory are reproduced and discuss two examples which illustrate 
the method and show possible advantages of this new way of constructing the Rayleigh- 
Schrodinger series. 

In section |^ the full equivalence of the two perturbation expansions in all orders is 
proven. 

Finally, in section ^ we discuss previous constructions mimicking classical perturbation 
expansions in quantum mechanics by Kummer [l^, Ali [Q, Eckhardt |T^, and Ben 



Lemlih and Ellison and their relation to the present method and conclude with some 
remarks about future investigations. 



2 Classical Poincare— von Zeipel perturbation theory 

In this section we will describe the classical Poincare-von Zeipel perturbation theory 
along with the method of averaging in a geometric manner such that its generalization 
to quantum mechanics is almost self-evident. To avoid later confusion we will use lower 
case letters for the classical situation. The unperturbed Hamiltonian ho is a function on 
phase space 7 which is equipped with a Poisson bracket structure {■,■}. ho is assumed 
to be sufficiently smooth, integrable and nondegenerate in the sense of Liouville- Arnold, 
i. e. it has n := ^diirt'^ functionally independent constants of motion bi, . . . ,bn which are 
in involution, define the invariant tori and have the property that 

{ho,g} = 0^ g = g{bi, . . .,bn). 
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Furthermore, let 



p=0 P- 



be the perturbed Hamiltonian where the perturbations hp, p > 1 are assumed to be 
sufficiently smooth functions on 7. The idea of Poincare-von Zeipel perturbation theory 
is to look for an e-dependent generating function 

00 

(with e-independent smooth functions wi) such that —w{e) generates a canonical flow 
(^(e) with "time" e. Then (p{e) := is a transformation on 7 determined uniquely by 

d 

— V?(e)* = ad w{e) o ip[e)* (3) 
^(0) = io?^ (4) 

and gives rise to the following action on phase space functions a: 

(p{e)*a := ao ip{e). (5) 

Here ad f{g) is defined for any two phase space functions f,gsis 

adf{g):={f,g} (6) 

and for future use we remark that {adw{e)y := adw{e) o ■ ■ ■ o adw{e) {p times). For 
later purpose we shall need an expansion of (p* in terms of some differential operators tp 
independent of e: 

^i^r = E z^tp- (7) 

p=o P- 

The tp are then recursively defined through to — id^ and 



'pi 



tp+i = E ( J adwi+i o tp-i. (8) 



Acting with (fi{e)* on the perturbed Hamiltonian h{e) gives a new Hamiltonian 

k{e) := ^(e)*/i(e) (9) 
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which is assumed to be analytic in e: 



OO 



and for which one finds 



He) = E ZlK (10) 



ko = ho (11) 
kp = adwp{ho) + fp, p>l (12) 

where /i = hi and for p > 2 

p-2 /p _ i\ 

fp-=hp + J2\^ ^ j {adwi+i{kp^i^i) + tp_/_i/i/+i). (13) 

If we had a solution iljk{t){t) : 7 — 7 of the motion with Hamiltonian k{e) then 

i^h(.)it) ■= ip{e) o V'fc(,)(t) o ip{e)-^ (14) 

would give us the desired solution for the perturbed Hamiltonian h{e). In general it is 
not possible to find Lp{e) such that a solution for k{e) may be found. However, we may 
choose the wi successively in such a manner that each kp is integrable, i. e. (since /iq is 
integrable and nondegenerate) such that 

adh(){kp) = Vp. (15) 

In this way one can trivially solve the equations of motion for k{e) up to any finite order 
in e and thus obtain via (|1^ solutions of h{e) up to the same order in e. 
Thus we have to find the Wp successively such that 

kp = adwp{ho) + fp (16) 
adho{kp) = 0. (17) 

Equation (|TB|) and (|T^ are solved by the method of averaging [|r^, |TB[. Let f3 = . . . , /3„) 
be the coordinates canonically conjugate to the integrals b = (61, . . . , 6„). Then the fiow 
ipho generated by the unperturbed Hamiltonian ho expressed in the coordinates (6, /?) is 



bo\^( Kt) \ / bo 

(3o i I Pit) ] I f3o + uj{bo)t 
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where 

gives the frequencies of the unperturbed motion which we assume to be independent 
over the rationals (non-resonant) for the given bo, i. e. c ■ uj{bQ) = ^ c = for any 
c = (ci, . . . ,c„) with integer Cj. Let g be any function on phase space which has the 
Fourier decomposition 

g{b,P)= J2 aibU""-^ (20) 

cSZ" 

and define the phase space functions 

g := hm - / (^,„(-t)*f7 (21) 

T^oo 1 Jo 

s{g) := lim - / dt ds (iphoisYg - g) (22) 

T^oo I Jo Jo 



then one finds 



and it is easy to see that 



9 = 9{h)o (23) 

s{g) = E (24) 



g = ad{s{g)){ho) + g (25) 
adhoig) = 0. (26) 



With this construction we choose now in (piq) for Wp'. 

Wp = siU) (27) 

then kp = fp commutes with ho as desired. The important point to note here is that we 
have formulated the basic (averaging) constructions (^1]) and (p2D necessary to solve ([T6| ) 
and {\[7\) in a coordinate free way. This geometric description using the time average is 
immediately suitable for generalization to self adjoint operators, i. e. to quantum me- 
chanics. It should be noted, however, that (^7|) is not the only possible solution of ([T6| ) 
and (0) since 

w'p:=Wp + Vp (28) 
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is also a solution of our problem as long as 

ad Vp{ho) = 0. (29) 

This nonuniqueness will also emerge in the quantum mechanical setting since it is also 
present in Kato's rigorous exposition of Rayleigh-Schrodinger perturbation theory [Q. 

3 Quantum Poincare— von Zeipel perturbation theory 
and averaging 

Now we shall develop the quantum mechanical analogue of the classical theory presented 
in the previous section. For this purpose we use capital latin letters to denote operators 
on some Hilbert space F. Let Hq be the unperturbed Hamiltonian operator which is 
assumed to be diagonalized in some basis and let 

oo p 

m := E -iH, (30) 

be the perturbed Hamiltonian. Here we do not impose any conditions (like e. g. bounded- 
ness) on the perturbations Hp, p > 1, but proceed purely on a formal level. A mathemat- 
ically rigorous justification of each step is notoriously intricate and will not be attempted 
here because it would completely obscure the basic ideas of the method. Since ultimately 
we shall prove the equivalence of this method to the usual Rayleigh-Schrodinger pertur- 
bation theory the conditions of the latter theory needed to guarantee convergence (see 
e. g. [|I^) will be sufficient to make the quantum Poincare-von Zeipel perturbation theory 
convergent as well. In analogy to the classical situation we seek a self adjoint generator 
(operator) 

oo p 

W{e) := E -l^P+i (31) 
p=o P- 

such that —W{e) induces the unitary flow H(e) with "time" e, i.e. 

^S(e) = ^W{e)E{e), 2(0) = 1. (32) 

Then $(e) := S(e)~^ is the unique solution of the initial value problem 

^$(e)* = ADW{e)o^ey (33) 
ae 

$(0) = 1 (34) 
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where $(e)* acts on any operator A via 

$(e)*A:= $(e)-M$(e) (35) 
and where ADF{G) is now defined as 

ADF{G):='-[F,G] (36) 



for any two operators F, G (again we omit the technical details necessary to make ( p6D well 
defined for unbounded operators) and as in the classical case {ADW{e)y := ADW{e) o 
• • • o ADW{e) {p times). 

As in the classical case it will be useful to expand $* in terms of e-independent 
operators Tp-. 

OO p 

$(6)* = Y: -T, (37) 
p=0 P- 

where the Tp are then recursively defined through Tq = 1 and 



Tp+i = E (yj Wi+i o Tp.i. (38) 



1=0 



Note that $(e)* and thus the Tp act on operators whereas $(e) itself is a transformation 
on Hilbert space which can also be expanded as 

OO p 

m = E -i^P (39) 
p=o P- 

and where the following recursive relation for the $p can be derived from (^21): 

Vi = -^E0V/W^^+i (40) 

and $0 = 1- 

Transforming the perturbed Hamiltonian H{e) with $(e) gives a new Hamiltonian 

K{e) := <l>{e)*H{e) = <l>{e)-^ H{e)<l>{e) (41) 
which is assumed to be analytic in e: 

OO p 

= E Zi^p (42) 
p=o P- 
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and for which one finds 



Ko = Ho (43) 
Kp = ADWp{Ho) + F„ p>l (44) 

with Fi = Hi and for p > 2 

p-2 



Fp := Hp + J2 I ^)(^^W^m(^p-^-i) + 7;-i-i^m. 



(45) 



All these equations are exactly analogous to the classical case but it is to be emphasized 
that they are perfectly well defined operator equations. But how are we to choose W 
now? If we could diagonalize K{e) up to a given finite order in e, we could read off its 
eigenvalues and eigenvectors to that order as well, but this means that we have found the 
eigenvalues and the eigenvectors of the perturbed Hamiltonian H{e) since by ( |^ ) H and 
K are unitarily equivalent. Before we write this out in formulae let us first see how we 
can diagonalize K order by order using the method of quantum averaging. It is obvious 
that the necessary and sufficient condition for diagonalization is the equivalent of (p!?!) , 
i. e. requiring 

ADHo[Kp) = Q, p>l (46) 

means that all Kp commute with Hq and thus Hq and Kp can be diagonalized simultane- 
ously, hence K can be made diagonal to any finite order in e. Consequently, in order to 
diagonalize K to any finite order we need to solve the quantum analogue of (|16|) and (p!7D, 
i. e. 

Kp = ADWp{ho) + Fp (47) 
ADHoiKp) = 0. (48) 

This, too, is done analogous to the classical case. Let be the unitary fiow generated 
by the unperturbed Hamiltonian Hq such that for any operator G 

^Hoi-ty G = exp(-^tiJo) G expi^Ho) =: G{t) (49) 



where the last equation introduces a simplified notation. Suppose now that G is such that 

G := lim - r dtG{t) (50) 

T^oo T Jo 

S{G) := hm 1 ^dt fds (G{s)--G) (51) 
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exist and such that 



then it follows that 



lim = (52) 



G = AD{S{G)){Ho) + G (53) 
ADHo(G) = 0. (54) 



We first prove (0): 



ADHoiG) = UHo,G]= hm UHo,Git)]dt 

tl T^QO I Jo tl 

= _ lini i r ^G{t)dt = - lim ^^^l' ^ = 

T^oo T Jo dt ^ ^ T^oo T 

by assumption. Moreover 

AD{S{G)){Ho) = lim- / dt ( ds UGis) ~G, Ho] = \im - [ dt f ds '-[G{s) Ho 

T^oo 1 Jo Jo n T-^oo 1 Jo Jo n 

hm - / dtf ds- G{s) = lim - / rft (G(t) - G) 

T->cxD _/ Jo Jo ds T-+cx> _/ Jo 



T- 

= G-G 

which proves (|53|). Equations (0)-(0) are the exact quantum analogue of the classical 
averaging technique with the noteworthy absence of any non-resonance condition. Thus 
for any p>l equation ( ^Tj) and (|^) are successively solved by 

W, = S{Fp) (55) 
Kp = Tp. (56) 



Using (1^) we may simplify F„ by noting that 



ADWi+^iFp^i^,) = ADWi+iiFp^i^,) (57) 



and with (|55|) we have Wi+i = S{Fi+i). Assuming continuity of the maps S and " one can 
formally show that S{B) = S{B) for any operator B for which S{B) and B exist. On the 
other hand it is evident that S{B) = 0. Putting these things together shows that 



ADWi+iiFp^i^^) = (5J 
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such that the expression ( |56D for Kp does not contain contributons arising from averaging 
the terms ADWi+i{Kp-i^i) in the expression (45) for the Fp. 

Let us summarize what we have done so far: Given a perturbed Hamiltonian H{e) = 
E^o have shown that chosing Wp = S{Fp) in W{e) = ^^p+i leads to 

K{e) = H, + Y.-Fp = $(e)-i H{e) ^{e) (59) 
p=i P- 

and the Fp all commute with the unperturbed operator Hq and are given by 

F\ = Hi. (60) 
Tp = TTp + Y^i^' jTp^i^,Hi+,, p>2. (61) 

The first few terms in the expansion (^9]) are 

Fo = 0, Ko = Ho, 

Fi = H,, Ki = H^ (62) 

F2 = H2 + ^[Wi,K, + H,], K2 = H2 + l[WuH,], 

^p=l p\ 

finite A^. Let us introduce the following notation 



where Wi = S{Hi). Hence, E^Li "-^Fp and Hq can be simultaneously diagonahzed for any 



N gp 



^'^(^) = T.-iFp (63) 
p=o Z^- 



which implies 

Let 
and 

Evidently 



^'^W = T.Zi% (64) 

p=0 

ir^(e) = (<l>^(e))-i i7(e) $^(e) + 0(e^+i). (65) 

^''(6)|jr(e) = ^f(e)|jr(e) (66) 

H{eme) = E,{e)\j){e) (67) 

i/(e)<l>^(6)|j)^(6) = i?f (e)$^(6)|j)^(e) + 0(e^+i), (68) 
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i. e. the eigenvalues Ej{e) and eigenvectors |j)(e) of the perturbed Hamiltonian H{e) are 
approximated as follows 

E,{e) = (6) + 0(6^+1) (69) 
b)(e) = <f"'(e)|jr(e) + 0(e^+^). (70) 

Hence, we have used the quantum analogue of the averaging method to construct a quan- 
tum mechanical perturbation theory. Just as in the classical case, however, the solutions 
for Wp constructed here are not the only ones. One encounters the same nonuniqueness 
as given by (pHf ) and (|2^) in the classical case. 



4 Examples: Discrete spectra 
4.1 General second order terms 

In this section we will apply the theory developed in section ^ to the case of a Hamiltonian 
Hq which is assumed to have a purely discrete spectrum with finite degeneracy 

Ho= \j,a)El{a,j\ (71) 

j; a&Dj 

where the sum over a runs over Dj := {l,...,dj = dim{Eig{Ho, Ej))}. For any self 
adjoint G one then obtains: 

G= Y \j,a){a,j\G\j,P){P,j\ (72) 

S{G) = - E \j^c.)^-^0^^{P,k\. (73) 

Using (§5), (0), and ( |73| ) one finds after straightforward calculations 

K\e) = E \j,a){E] + e{a,j\H,\j,P)){(3,j\ (74) 

+ E M ^ + L F0_ FO 

j; a,l3eD, \ [ ^ j^k; ^dD^ ^3 ^k } J 

Consequently the eigenvalues Ej^^{e) of K^{e) are determined as solutions of the secular 
equation of dj- dimensional matrices: 

det [E^6ap - {a,j\K\e)\j,P)) = (75) 
12 



which coincides with the usual Rayleigh-Schrodinger result. The corresponding eigenvec- 
tors of -^^(e) are then 

U,«)'(e)= E (76) 
and one has for the eigenvector |j, a)(e) of H{e) 

b, «) (e) = (l - + ^ ((^^i)' - ^^2) ) I J, a)\e) + 0(6^) (77) 

with 

Wi = - 2^ — -g (78) 



and 14^2 = '5'(i/2 + ^[W^i; Hi + //i]) which we shall not write down here but which inserted 
into (|7^) yields the corrections to the eigenvctors to second order known from Rayleigh- 
Schrodinger perturbation theory in the case of non-degenerate spectrum. In this case we 
see from (76) that the eigenvectors of i^^(e) and Hq coincide: |j)^(e) = In fact, since 



K^{e) = K'^-\e) + —Kr, (79) 

and [Hq, K]^] = VA^ we have by induction 

|j)^(6) = \j) VAT (80) 

in the non-degenerate case. 

The general equivalence of the quantum version of the Poincare-von Zeipel perturba- 
tion theory to the standard Rayleigh-Schrodinger perturbation theory will be proven in 
section ^ The fact that the corrections to the eigenvalues in the Poincare-von Zeipel the- 
ory are derived from an averaging procedure may, however, provide some computational 
advantage since it gives the sums over intermediate states so common to standard pertur- 
bation theory in closed form. For example, the second order term in the non-degenerate 
case is given by 
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4.2 Example 1: Anharmonic oscillator 

For the nondegenerate case we shall illustrate the method in the example of the harmonic 
oscillator Hq = ^ + with cubic perturbation Hi = jx^, Hp = 0, p > 2 (anhar- 

monic oscillator with ^ = 1) where the quantum Poincare-von Zeipel method will permit 
us to compute the sums over intermediate states and the corrections up to O(e^) without 
much effort. All calculations are straightforward if we use the operators 

^/2^dx^^^ '~ a/2 ( ^^^^ 

for which one finds 

a{t) = e**a a\t) = e-''o) (83) 

such that 

= ^Ae'''a' + e-''\a^Y + Ae'''aH,a + Ae-'''a'^H,a^) + l{H,f + ^ 

Wx = -^((a^Y -a'^ + Sa^Hoa^ -SaHoo) 
64z ^ ^ 

Wi{t) = ( e-^'*(a^)" - e^'*a^ + 8e~^''a^Hoa^ - Se^^'aHoo) 

K2 = ^{[{a^)\a''] + 32[a^Hoa\aHoa]) 

^ .1 /3, .2 3 \ 2 .3 51 .2 59 . 21 



' '2 VS'-" " 16/ V64'' 12S' 12S' 128 

This formula for the correction to the eigenvalues was first derived by Heisenberg and 
is also reproduced by Kummer using his normal form approach ||12|| . 

4.3 Example 2: Henon— Heiles system 

To illustrate the method for the degenerate case we apply it to the two-dimensional 
Henon-Heiles system whose unperturbed Hamiltonian is (we choose again ^ = 1) 



I ( £ (P\ (Xi)2 + (X2f 1 1 
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where the Nj, j = 1,2 are the number operators and the perturbation is 

Hi{a, P) := a{xifx2 + /3(x2)^ Hp^O if p>2. (85) 

It is only for convenience that we have chosen the "degenerate" case (i. e. equal frequencies 
for the two one-dimensional oscillators in Hq). The method is completely oblivious to 
that distinction. In this example we treat a and /? as one perturbation parameter in the 
sense that a = ea, j3 = e(3 and e is the single perturbation parameter which we set equal 
to one at the end of the calculation. As in the case of the anharmonic oscillator it is very 
convenient to use the operators 

which also evolve according to 

aj{t) = e^*aj a]{t) = e-'^a]. (87) 

This yields 

Hi{t) = -l={e3^*(a(ai)V + /3(a2)^)+e-=^^*(a(aI)M + /?(4)') 

+ e'' (2aHoia2 + a(ai)^4 + 2/3iJo2a2 + /9(a2)^4) (88) 

Ih ^ Q (89) 

and 

W^ = ^{^(a(al)M + /3(4)V^H«i)'«^ + '^M') 

+ q;(4)^ + 2Q;i/oi4 + /3(4)^a2 + 2/3//024 (90) 
- {2aHma2 + a{aifal + 2pHo2a2 + /3(a2)^4) } 

From this one can read off Wi{t) and obtains after some tedious but straightforward 
calculations (which we have executed with the help of the symbolic computation language 
MAPLE) 

^2(«,/3) = -^-^ + Y^ + -^ + (ai)'(4)' + (a2)Vi)' + ^ + ^ 
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2 2 



(91) 



+ 3 A^2 + 3 A^i a/?. 



Keeping in mind that we set e = 1 and that = i^i = we now have to find the 
eigenvalues E"^^ f^) -^^(«, /3) = -ffo + ^-^2 which will give us the correct eigenvalues 
of Hq + Hi up to second order. Let k = denote the ground state (no degeneracy: (0, 1)), 
k = 1 the first exited state (double degeneracy: (1, 1), (1, 2)), and k = 3 the second exited 
state (triple degeneracy: (2, 1), (2, ±)) then we find the following eigenvalues 



^?o,i)(«,/5 



1 
2 
2 
3 
3 
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3 /5q! 


8 


24 


4 


11/3^ 


11^2 


9/?q; 


8 


8 


4 


71 /?2 


13^2 




8 


24 


4 


71/32 


19 


27/3a 


8 


8 


4 


101/32 


15/3 


a 17 



(92) 



± 



4 8 4 

The results for E"^ agree with those obtained by Kummer and Ali (except for the 

^(2,1) 



factor of q;2 in i??2 d in which is probably due to a typographical error) 



5 Equivalence to Rayleigh— Schrodinger perturbation 
theory 

5.1 Non— degenerate case 

In section ^ we have already seen that at least up to second order the quantum analogue 
of the Poincare-von Zeipel method and the standard Rayleigh-Schrodinger perturbation 
theory coincide. In this section we show that this is indeed true for the full perturbation 
expansions. To do this we recall briefiy how the standard Rayleigh-Schrodinger expansion 
is constructed in the nondegenerate case. With the help of a suitably chosen contour 
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integral in the complex i?-plane one can show that the projector 

1 r dE 



on the jth eigenspace of H{e) is analytic in e and that for e sufficiently small {j\Pj{e)\j) > 
0, 0, |l^. This projector then gives a normalized eigenvector |j)(e) of H{e) to the eigen- 
value Ej{e) via 

From this one obtains 

^ mmmi^) ^ <^i^|;}g;:>;;f"i^> . (95) 

Using the expansion for Pj{e) the right hand side of (p5|) then yields an expansion for 
Ej{e) in e which is the usual Rayleigh-Schrodinger perturbation series. 
On the other hand it follows from (^) that 

Kie)\j) = E,{e)\j) (96) 

where K{e) = $(e)~^ H{e) $(e) which implies 

m^) = m\j) (97) 

Inserting ( p8[ ) in ( ]95| ) yields 

oo p 

= m{e)\j) = T.-mAj) (99) 

p=o ^'• 

which proves the equivalence in all orders for the non-degenerate case. 
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5.2 Degenerate case 

Suppose now that HqU, a) = Ej\j, a) with possible degeneracies a G Dj = {1 , . . . , dj} , 
that Ej is an isolated point of the spectrum (t{Hq), that Pj = Y^a^Dj \ is the 

projector on the jth eigenspace of Hq and that Pj{e) defined as in ( P5| ) exists and is 
analytic in e. Then it has been shown that @, [T^ 



[Hie) \RanP,(e)) = ^(^(e)) {E \ \E ~ E^\ < r} (100) 

(where A luanB means restriction of the operator A to the range of B) and that there 
exists a unitary operator U (e) such that 



and 

satisfies 



P,(e) = f/(e)P,^(e)-^ (101) 

^(6) := U{e)-'H{e)U{e) (102) 

H{e)P,=P,H{e). (103) 

Then Pj{H{e) — E)Pj is a finite dimensional operator analytic in e and the eigenvalues 
Ej^ai^) of H{e) are found as the dj roots of the equation 

det{P,{H{e)-E)P,)=0. (104) 

The operators Lj(e) := PjU{e)~^ and Rj{e) := U{e)Pj sandwiching H(e) in ( |104| ) satisfy 
certain differential equations involving Pj (e) whose expansion in terms of e is known from 
(P^). With the help of these differential equations Lj{e) and Rj{e) can be expanded in e as 
well. The A^th order approximation in the Rayleigh-Schrodinger series for the degenerate 
case is then obtained by solving ( |104| ) where terms of order higher than N are neglected. 



We refer the reader to [g for more details on Kato's rigorous exposition of the usual 
quantum mechanical Rayleigh-Schrodinger perturbation theory. 

As pointed out in the unitary transformation [/(e) exists but is not necessarily 
unique. As we shall see below this nonuniqueness is equivalent to the nonuniqueness of 



the choice of Wp mentioned at the end of section 



In order to prove the equivalence with the quantum Poincare-von Zeipel method we 
first note that 

m-'U{e)P,=P,m-'U{e). (105) 
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In fact, since [K{e),Ho] = we may write the orthonormalized eigenvectors of K{e) as 
J2/3eDj '^''{^)ai3\j, P) where u^{e) is a (i-,— dimensional unitary matrix. Hence, the projector 
on the space ®a&D Eig{K{e), Ej^ai^)) is 



Y: u^ieU\j,P){7,j\{u^ie))\ = Y: \3,P){P,3\ = Pr (106) 

Since K(e) = $(e)"^i/(e)$(e), one has 

P,(e) = <l>(e)P,$(e)-i (107) 
which together with ( |101| ) proves ( |105| ) . Equation ( |105| ) states that 

Z(e) := $(e)-i f/(e) (108) 
is a unitary transformation commuting with all Pj. On the other hand, from 

$(e) K{e) $(e)"^ = H{e) = U{e) H{e) U{e)-^ (109) 

it follows that 

H{e) = Zier'Kie)Zie) (110) 
which implies that the roots of equation ( |104D are identical to the roots of 

detiP,iKie)-E)P,)=0 (111) 

and this proves the equivalence of the two methods for the eigenvalue-expansions in the 
presence of degenerate eigenvalues. 

Moreover, as can be seen from (|110|) the eigenvectors of H{e) are related to those of 
K{e) by the unitary transformation Z{e) which preserves each eigenspace of Hq. It is 
likely that the nonuniqueness in the choice of Wp may be exploited to make Z{e) trivial 



11] but this is still under investigation. 



6 Discussion and conclusion 



Kummer |12] was the first to discuss the averaging method for quantum systems. Based 
on ideas from classical averaging he constructed a perturbation method called the normal 
form approach |jl3| which is equivalent to "time averaging" but instead of using 
averaging to solve (jS7\) and (H3) it employs algebraic constructions. 
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Motivated by the Birkhoff-Gustavson normal form in classical mechanics Ali []T^ has 
developed a quantum analogue of this and his construction yields the same expansion as 
that of Kummer. Working explicitly with an algebra of destruction and creation operators 
Eckhardt [1^ has also constructed a quantum analogue of the Birkhoff-Gustavson normal 
form. A quantization of the classical Birkhoff-Gustavson normal form has been attempted 
by Robnik p^] but this is necessarily plagued by ordering problems which do not affect 
our work and the other contributions cited above (this is only partly true for [ll5|]). 

The constructions of Kummer, Ali and Eckhardt have in common that the existence 
of the generators of the unitary transformation has to be assumed or assured by certain 
additional conditions. In the present paper the necessary generators Wp are (at least 
formally) explicitly constructed. 



In fact, it can be shown |TI| that the method of quantum averaging as presented 
here provides explicit solutions in terms of the time averaging integrals for the algebraic 
constructions of Kummer. The algebraic constructions have the advantage of rigorous 
validity but lack constructive procedures needed to execute the algorithm. The approach 
presented here has - apart from its conceptual proximity to the classical situation - the 
advantage of providing explicit constructions. Due to the analytic character of these 
constructions, however, technical problems which are absent in the algebraic approach 
may arise. 

Using a slightly modified quantum averaging in our sense for the particular example 
of the quantum anharmonic oscillator Ben Lemlih and Ellison have derived rigorous 
error bounds on approximations to the quantum time evolution. 

Their work also contains a suggestion to compare the approximation to the eigenvalues 
of this specific problem to the usual perturbative corrections, i. e. Rayleigh-Schrodinger 
theory. Ali and and Kummer |]l3l have found that in all the examples they have 
treated the normal form perturbative results agree with Rayleigh-Schrodinger theory 
(incidentally this is not true for Robnik's expansions |2^) and Eckhardt also suggests 
that the Birkhoff-Gustavson perturbation expansion is identical to the usual Rayleigh- 
Schrodinger pertubation theory. Our work then provides an explicit proof of this assertion 
since, as Kummer has shown, his normal form approach is equivalent to the averaging 
method in quantum mechanics and we have shown that averaging is completely equivalent 
to the Rayleigh-Schrodinger theory yielding the sums over intermediate states in closed 
forms. 

A very important aspect related to this work (but not discussed here) is the fact, that 
just as in classical mechanics a superconvergent perturbation theory can be constructed 
with the help of averaging, this can be done in quantum mechanics as well and yields a 
perturbation theory explicitly distinct from the usual Rayleigh-Schrodinger theory H, 
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Work is in progress to establish the technical conditions necessary to put quantum 
averaging on a rigorous mathematical footing and to determine how the nonuniqueness 
of the Wp may be used to trivialize Z{e). It may also be possible that standard time- 
dependent perturbation techniques (e. g. sudden aproximation) can be formulated as 
analogues of classical time-dependent averaging. 
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